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We present a method to measure the weak phase 2(3 + 7 in the three-body decay of neutral B°(B°) mesons to the final 
states D K°^- These decays are mediated by interfering amplitudes which are color-allowed and hence relatively large. 
As a result, large CP violation effects can be observed with high statistical significance. In addition, the three-body 
decay helps to resolve discrete ambiguities present in measurements of the weak phase. The experimental implications of 
conducting these measurements are discussed, and the sensitivity of the method is evaluated using a simulation. 



r- ■ 
m ■ 

: 

o . 
m ■ 

o ; 

^ : 

a: 
1 

<D . 
■ 



X 



1 Introduction 

The violation of the CP symmetry is now established 
in the B meson sector and the parameter sin(2/3) is 
measured with precision by BaBar [1] and Belle [2]. 
However, the measurement of the other angles of the 
unitarity triangle are necessary for a more compre- 
hensive study of CP violation. These measurements 
suffer from several difficulties related to the uncer- 
tainties of the theory, the limited statistics and/or 
ambibuities in their extraction. We have recently pro- 
posed a new method [3] involving three-body B meson 
decays, which circumvent these problems and allows to 
determine the angle 7 = a,rg(-V ud V* b /V cd V* b ). 

In this paper we extend this method to the weak phase 
Mg(V t fV t 2 b V: b V us /V c * d 2 V cb V cs ), which can be identi- 
ficd to the combination 2/3 + 7 of the unitarity triangle 
angles [4] . Important constraints on the theory will be 
obtained from the measurement of this phase. In the 
following, we describe a new time dependent method 
involving the entire area of the Dalitz plot of the 3- 
body decays B° -> D ± K^. 

The method involves no reconstruction into CP eigen- 
states and avoids problems with interfering DCSD. 



2 Color-Allowed B°^D ± K^ decays 

We investigate a way to circumvent the color suppres- 
sion penalty by using B°(B°) decay modes which po- 
tentially offers significantly large branching fractions 
and CP asymmetries. The particular decays consid- 
ered here are of the type B° -> D { *^ tt / p ± . 
These three body decays may be obtained by popping 
a qq pair in color allowed decays. Although modes 
where one or more of the three final state particles is 
a vector can also be used, for clarity and simplicity 
only the mode B° — > D^K^n^ is discussed here. 




Figure 1. Feynman Diagrams for the decay B° — > 
D ± K^tt^ 1 involving the CKM matrix element product 

v; b v us or v: b v ca . 



The diagrams leading to the final states of interest 
are shown in Fig. 1. As can be seen, the diagrams 
(Fig. la and lc) are both color-allowed and of order 
A 3 = sin 3 (0 c ) in the Wolfenstein parameterization [5], 
where 9 C is the Cabibbo mixing angle. 

The general formalism describing how the weak angle 
is extracted from observation of the time-dependent 
CP violation has been first presented in [6]. How- 
ever, in the present analysis additional information is 
used to resolve the inherent ambiguities related to the 
method. 

Selecting a particular point i in the Dalitz plot, let us 
write the amplitudes corresponding to the transitions 
in Fig. la and b on the one hand and lc and d on the 
other hand as 



Ai(B° -» D-K° s tt+) 
Ai{B° -> D+K°w-) 



A Cl e lSc 



(1) 

(2) 



2 
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and their CP conjugates 

Ai(B° -» D+K° s ir-) = 
Ai(B° -> D~ K° s -k+) = 



A Ci e 



(3) 
(4) 



where 7 is the relative phase of the CKM matrix ele- 
ments involved in this decay, and Ac (An) and 8 c (&u) 
are the real amplitudes and CP-conserving strong in- 
teraction phases. 

The time dependent probability density function 



(PDF) of the transition B°~> D-K s n+ i 



is 



t (-) 



(-) 
1 + 



IZi cos x— + T>i sin (2/3 + 7 + ASA sin x— 

T T 



(5) 



where x = Amr, 7£j 



2 Pi 



with pi = and A£j = <5ci — <5wi- We have used the 

approximation \p/q\ = 1. ( 

The PDF of the transition E°— > D+ K%~ is 



(-), 



I B ^ D + K%- I = 



1 



(+) 



IZi cos x— + sin (2/3 + 7 — AS{) sin x 
r 



(6) 



From the total number of events in each of the four 
final states (5) and (6) and a global fit of their time 
dependence, it is possible to extract the quantities, 
A c 2 i + A u \ , ft , S , S where S = sin(2/3 + 7 + A<5 t ) and 
S = sin(2/3 + 7 - ASi). One obtains: 

sin 2 (2/3 + 7) = \ [\ - SS ± ^ 2 )(1 - S*)) (7) 

Hence, in the limit of very high statistics, one would 
extract sin 2 (2/3 + 7) for each point i of the Dalitz 
plot. However, for every point of the Dalitz plot, 
2/3 + 7 is obtained with an eight-fold ambiguity in 
the range [0, 2n], because of the invariance of the 
sin(2/3 + 7 ± ASi) terms in Eq. (5) and (6) under the 
three symmetry operations 



Sir/2 



2/3+7 -> 7r/ 2 - M As 
2/3+7^ -(2/3+7), A«J 
2/3+7 -» 2/3+7 + t, A<5 



7T/2 - (2/3 + 7) 

7T-A5 (8) 

A5 + 7T 



When the multiple measurements made in different 
points of the Dalitz plot are combined, some of the 
ambiguity will be resolved, in the likely case that the 
strong phase ASi varies from one region of the Dalitz 
plot to the other. This variation can either be due 
to the presence of resonances or because of a varying 



phase in the non-resonant contribution. In this case, 
the exchange symmetry S n /2 is numerically different 
from one point to the other, which in effect breaks this 
symmetry and resolves the ambiguity. Similarly, the 
S- symmetry is broken if there exists some a-priori 
knowledge of the dependence of ASi on the Dalitz plot 
parameters. This knowledge is provided by the exis- 
tence of broad resonances, whose Breit-Wigner phase 
variation is known and may be assumed to dominate 
the phase variation over the width of the resonance. 
To illustrate this, let i and j be two points in the 
Dalitz plot, corresponding to different values of invari- 
ant mass of a particular resonance. One then measures 
sin(2/3+7±A5j) at point i and sin(2/3+7±(A5 i +a i j)) 
at point j, where is known from the parameters of 
the resonance. It is important to note that the sign 
of a>ij is also known, hence it does not change under 
S-. Therefore, should one choose the S- -related so- 
lution sin(7r — (2/3 + 7) + ASi) at point i, one would 
get sin(7r — (2/3 + 7) + (ASi — a ij)) & t point j. Since 
this is different from sin(2/3 + 7 ± (ASi + a ij)) the S- 
ambiguity is resolved. This is illustrated graphically 



in Eq. 9, where = 2/3 + 
sin(</> ± ASi) ^ 
BW I 

S- 

sin(0± (ASi+aij)) ^-/-^ 



7 is the weak phase: 
sin(7r — (f> + ASi) 

I BW (9) 
sin(7r - 4> + (ASi — aij)) 



Thus, broad resonances reduce the initial eight-fold 
ambiguity to the two-fold ambiguity of the S„ sym- 
metry, which is never broken. Fortunately, S„ leads 
to the well-separated solutions 2/3 + 7 an d 2/3 + 7 + n, 
the correct one of which is easily identified when this 
measurement is combined with other measurements of 
the unitarity triangle. 

3 The Finite Statistics Case 

Since experimental data sets will be finite, extracting 
7 will require making use of a limited set of parame- 
ters to describe the variation of amplitudes and strong 
phases over the Dalitz plot. The consistency of this 
approach can be verified by comparing the results ob- 
tained from fits in a few different regions of the Dalitz 
plot, and the systematic error due to the choice of the 
parameterization of the data may be obtained by using 
different parameterizations. 

A fairly general parameterization assumes the exis- 
tence of Nr Breit-Wigner resonances, as well as a non- 
resonant contribution: 



At(B°^D-K°n+) 



(a Cq e w <=° + Acj B Sj (£) e^) 



,iSc(£) 



A ( (B° -» D+Ky-) 



(10) 



,*(M£)+7) 
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where £ represents the Dalitz plot variables, B Sj (£) = 
b Sj (£) e lSj ® is the Breit-Wigner amplitude for a parti- 
cle of spin Sj, normalized such that J \b a A£)\ 2 d£ = 1, 
and (5^o (^co an d ^co) are the magnitude and 
CP-conserving phase of the non-resonant b — ► mcs 
(6 — > cus) amplitude, and ,<4#j and 5uj (Acj and <5cj) 
are the magnitudes and CP-conserving phase of the 
b — » Tics (6 — > cms) amplitude associated with reso- 
nance j. The functions <5c(£) an d <fo(£) mav be as- 
sumed to vary slowly over the Dalitz plot, allowing 
their description in terms of a small number of pa- 
rameters. The decay amplitudes of B° mesons are 
identical to those of Eqs. (10), with 7 replaced by —7. 

The decay amplitudes of Eqs. (10) can be used to con- 
duct the full data analysis. Given a sample of iV e 
signal events, 7 and the other unknown parameters of 
Eq. (10) are determined by minimizing the negative 
log likelihood function 

N e 

X 2 = -2glogP(&), with P(0 = \Az(f)\\ (11) 
i=i 

where the amplitude A^(f) is given by one of the ex- 
pressions of Eqs. (10), or their CP-conjugates, depend- 
ing on the final state /, and are the Dalitz plot 
variables of event i. 

In what follows, we discuss important properties of the 
method by considering the illustrative case, in which 
the b — ► Tics decay proceeds only via a non-resonant 
amplitude, and the b — > cus decay has a non-resonant 
contribution and a single resonant amplitude. For con- 
creteness, the resonance is taken to be the _?T* ± (892). 
We take the a priori ^-dependent non-resonant phases 

to be a c (0 = MO = o. 

4 Simulation Studies 

To study the feasibility of the analysis using Eq. (11), 
we conducted a simulation of the decays B° — > 
D^K^ir^ and B° — > D ± K°ir T . Events were gener- 
ated according to the PDF in Eqs. (5) and (6), with 
the parameter values given in Table 1. In this table 
and throughout the rest of the paper, we use a tilde to 
denote the "true" parameter values used to generate 
events, while the corresponding plain symbols repre- 
sent the "trial" parameters used to calculate the ex- 
perimental x 2 - F° r simplicity, additional resonances 
were not included in this demonstration. However, 
broad resonances that are observed in the data should 
be included in the actual data analysis. 

The simulations were conducted with a benchmark in- 
tegrated luminosity of 400 fb _1 . The reconstruction 
efficiencies were based on current T(4S) detector ca- 
pabilities. We assumed an efficiency of 50% for recon- 
structing the K a , and 90% for reconstructing the 7r ± . 



Parameter 


Value 


Parameter 


Value 







Auo/Ago 


0.4 


MO 





Ack'/Aco 


1.0 


5k. 


1.8 


Ack* 


V2 x io- 4 r B 


Sua 


1.0 


2/3+7 


2.00 



Table 1. Parameters used in the simulation. The value of 
Ack* is chosen so as to roughly agree with the measure- 
ment of the corresponding branching fraction [7]. 



The product of reconstruction efficiencies and branch- 
ing fractions of the D^, summed over the final states 
if :F 7r ± 7r ± and is taken to be 4%. Finally the 

analysis efficiency including perfect tagging was esti- 
mated to be 20%. The numbers of signal events ob- 
tained in each of the final states with the above ef- 
ficiencies and the parameters of Table 1 are listed in 
Table 2. 

In Fig. 2 we show the \ 2 dependence as a function 
of 2/3 + 7 an d Sua- The smallest value of y 2 is 
shown as zero (white). At each point in these fig- 
ures, the x 2 is calculated with the generated values 
of the amplitude ratios Au /Acq = Auq/Aqq and 
Ack* I Acq = Ack* I Acq- When these amplitude ra- 
tios are determined by a fit simultaneously with the 
phases, the correlations between the amplitudes and 
the phases are found to be less than 10%. Therefore, 
the results obtained with the amplitudes fixed to their 
true values are sufficiently realistic for the purpose of 
this demonstration. 




It/2 it 3it/2 2>t 



Weak phase 2p+y 

Figure 2. \ 2 as a function of 2/3 + 7 and 5u , with the 
parameters of Table 1. 

In Fig. 3 we show the one-dimensional minimum pro- 
jection x 2 ( 2 /3 + 7) = min{x 2 (2/3 + 7, 6 Uo )}, showing 
the smallest value of x 2 for each value of 2/3 + 7. 



4 
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1D Likelihood scan 




n/2 n 3x12 2n 

Weak phase 2p+y 

Figure 3. Minimum projection of x onto 2/3 + 7. 



With equal resonant and non-resonant b — > cus am- 
plitudes, the ambiguities originating from the non- 
resonant regime are dominant, due to the relative sup- 
pression of the resonant interference terms, as they 
overlap only little with the non-resonant amplitudes. 
These are however significant enough to resolve all but 
the SV ambiguity. 
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Figure 4. The error on 2/3 + 7, 0-2/3+-, 
(top) 2/3+7 an d (bottom) Auo/Aco- 



as a function of 



In Fig. 4 we present 02/3+7, the statistical error on 
2/3 + 7 as a function of 2/3+7 an d -Au /Acq, obtained 
by fitting simulated event samples. Each point is ob- 
tained by repeating the simulation 250 times, to min- 
imize sample-to-sample statistical fluctuations. All 
the parameters of Table 1 were determined by the 



Mode 




N 

1 v ev 


Mode 






N ev 




D~K^TT + 


112 


B^ 


£H 




33 






111 


~B° — > 


D~ 




33 



Table 2. The numbers of perfectly tagged events of each 
type expected in 400 fb _1 obtained by averaging 100 sim- 
ulations using the parameters of Table 1 and the recon- 
struction efficiencies listed in the text. 

fit, and the arrow in these figures indicates the point 
corresponding to the parameters of Table 1. One ob- 
serves that 02/3+7 does not depend strongly on the 
value 2/3 + 7. As expected, a strong dependence on 
Auq/Aco is seen. However, it should be noted that 
02/3+7 changes very little for all values of Auq/Aqq 
above 0.4. This suggests that the potential for a sig- 
nificantly sensitive measurement is high over a broad 
range of parameters. 

5 Conclusions 

We have shown how 2/3 + 7 may be measured in the 
color-allowed decays B° — > D^*'^ K^*' 0, k / ' , focusing 
on the simplest mode B° — > D^K^it^. The absence 
of color suppression in the b — > ucs amplitudes is ex- 
pected to result in relatively large rates and signifi- 
cant CP violation effects, and hence favorable exper- 
imental sensitivities. While the Dalitz plot analysis 
implies some experimental complication, it should re- 
duce the eight-fold ambiguities, which are a serious 
problem with other methods for measuring 2/3 + 7. As 
a result of these advantages, this method is likely to 
lead to a measurement of 2/3 + 7 even with the current 
generation of B factory experiments. 
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